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1.  INTRODUCTION 


The  dynamic  response  of  buried  structures  subjected  to  explosive  loadings 
can  be  analyzed  using  many  different  models.  These  range  from  simple  models 
that  involve  a  few  parameters  to  extremely  complicated  models  that  represent 
many  details  of  the  explosion  process,  the  surrounding  ground  materials,  and 
the  structure  itself.  The  difficulties  and  uncertainties  associated  with 
simultaneously  analyzing  the  explosion,  the  ground,  and  the  structure  usually 
lead  to  consideration  of  procedures  that  allow  detailed  study  of  some  parts 
of  the  problem  while  treating  the  rest  in  a  simplified  fashion.  Approximate 
procedures  of  this  type  will  be  referred  to  as  uncoupling  methods. 

The  main  purpose  of  this  study  is  to  consider  an  uncoupling  method  that 
permits  detailed  analysis  of  the  structure  while  the  ground,  which  may  exhibit 
hysteretic  behavior,  is  represented  by  a  relatively  simple  model  based  on 
an  extension  of  the  plane  wave  approximation.  This  method  is  a  special  case 
of  a  more  general  approach  (currently  under  development  by  Weidlinger 
Associates,  and  is  summarized  in  the  Appendix)  to  uncouple  the  motion  of 
a  structure  from  that  of  a  nonlinear  medium  with  which  it  is  in  contact. 

The  plane  wave  approximation  has  been  used  as  an  uncoupling  method  for 
linear  fluid-structure  interaction  problems  for  many  years  (e.g.  see  Ref.  [2]). 
Recently,  Ref.  [3],  it  has  been  applied  to  the  nonlinear  problem  of  c  structure 
submerged  in  a  cavitating  fluid.  It  has  also  played  a  role  -  In  the  form  of 
viscous  dampers  -  in  the  development  of  "transmitting"  boundaries  for  soil- 
structure  interaction  calculations  (both  linear  and  nonlinear);  here  however 
the  boundary  is  placed  relatively  far  from  the  structure.  Jts  application  as 
a  soil-structure  uncoupling  method,  wherein  the  soil  is  entirelv  replaced  bv 
the  plane  wave  approximation  at  the  soil-structure  interface,  appears  to  be  quite 
recent.  It  was  successfully  applied  in  this  way  to  a  linear  two-dimensional 

problem  in  Ref .  [4]  and  to  a  nonlinear  two-dimensional  situation  in  Ref.  [  ">  ] . 
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2.  ONE  DIMENSIONAL  SOIL-STRUCTURE  INTERACTION 
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Consider  the  soil-structure  interaction  problem  illustrated  in  Fig.  1. 

A  possible  one  dimensional  model  of  the  roof  of  the  structure  in  contact  with  the 
soil  is  shown  in  Fig.  2.  This  one  dimensional  problem  will  be  taken  as  the 
model  with  which  the  uncoupling  method  will  be  compared  in  detail. 

2.1  Formulation  of  the  Continuum  Equations  (see  Fig.  2) 

The  soil,  having  density  and  depth  H^,  is  represented  by  piecewise 
linear  stress-strain  behavior  with  initial  loading  modulus  My,  and  unload¬ 
reloading  modulus  My,  where  My  ^  My  .  The  structure,  having  density 

and  depth  H^,  is  represented  by  a  linear  elastic  stress-strain  behavior 


with  modulus  M 


2  * 


The  shearing  resistance  of  the  roof  is  represented  by  an 


elastic,  ideally  plastic  continuous  support  that  exerts  a  force  F<,  where  the 
elastic  force-displacement  behavior  is  represented  by  a  spring  constant  K  and 
the  plastic  limit  is  represented  by  a  force  Fq  . 

Let  x  denote  a  coordinate,  increasing  with  depth,  with  x  =  0  corresponding 

to  the  ground  surface.  Let  u(x,t)  be  the  particle  displacement  at  point  x 

3u 

at  time  t  in  the  continua,  and  v(x,  t)  =  be  the  corresponding  particle 
velocity.  The  stress,  negative  in  compression,  is  denoted  by  a(x,  t).  Ignoring 
convective  effects,  the  continuum  equations  for  the  soil  can  be  written  as 


By. 

m 


.  3v  3  a 

pi  h‘5 


3a  _  3v 
3t  1  3x 


(1) 

(2) 


I 


where  My  is  either  My  or  My,  depending  on  the  state  of  the  material.  The 
continuum  equations  for  the  structure  can  be  written  as 


l*; 

fe 

Vs 


.  3v.3fl.  s 

2  3t  3x  l^A 
3a  3v 

-  =3  M  — — 

3t  2  3x 


(3) 

(4) 


where  F  depends  on  the  state  of  the  shear  support,  and  A  is  the  area  of 

w 

the  structure  that  is  in  contact  with  the  soil. 


2.2  Finite  Difference  Equations 


Equations  1  to  4  are  solved  numerically  by  finite  differences  using 

central  differences  in  space  and  time  on  a  staggered  grid  of  stress  and 

velocity  points,  as  shown  in  Fig.  3.  Stress  points  are  placed  at  x  =  0  and 

x  =  (H^  +  i^)  and  a  velocity  point  is  placed  at  x  =  .  In  the  soil  region 

difference  equations  are  expressed  in  terms  of  the  lengtn  Ax^  =  H^/(N^  +  1/2) 

where  is  an  arbitrary  positive  integer.  Similarly  in  the  structure  the 

difference  equations  are  expressed  in  terms  of  the  length  Ax2  =  ^/(^  +  1/2) 

Let  0^  denote  o(x t)  where  x^  =  (i-l)Ax^  for  i  =  1,  2,...,  +  1  and 

x.  =  H  +  (i  -  N  -  3/2) Ax„  for  i  =  N.  +2,...,  N  +  N  +  2.  Let  v  denote 
1  1  1  2  1  Li  l 

v(xi>  t)  where  x±  =  (i  -  l/2)Ax1  for  i  =  1,  2,...,  N1  +  1  and  xi  =  H1  + 

(i  -  Nj.  -  l)Ax2  for  i  «  +  2,...,  NL  +  N2  +  1#  Difference  approximations 

to  Eqs .  1  and  2  can  then  be  written  as 

vn+1  =  vn  +  -  (0"  -  0?)  for  i  =  1,  2,...,  N  (5: 

1  1  PjAx^  1+i  1  1 

M,At 

n+1  n  ,  1  ,  n+1  n+1.  ...  ,  .  / r  ■ 

0  =  0.  +  -7 -  (v.  -  v.  ,)  for  1  =  2,...,  N-.  +  1  (o, 

i  1  Ax^  1  i-l  1 

where  the  superscripts  denote  time  tn  or  t^+^  and  At  =  tn+1  -  t^  is  the 


time  increment . 


The  symbol  in  Eq.  5  denotes  a  spatially  averaged  particle  velocity, 


defined  by 


nn.  ,  n  n,  .  n  n  , 

V.  =  V.  +  CX  (v-,i  ~  v-)  “  a,-  i/o  (v.  -  v.  ,) 

1  1  1+1/2  l+l  l  1-1/2  i  l-l 


where  every  a  must  be  non-negative.  Small  positive  values  of  a  introduce 
a  linear  (artificial)  viscosity  which  can  be  used  to  suppress  spurious 
high-frequency  oscillations  that  are  produced  by  this  numerical  method. 

Similarly,  difference  approximations  to  Eqs.  3  and  4  can  be  written 


n+l  ~n  At  ,  n  n  ,  At  n 

Vi  Vi  P2&X2  *  P2H2A  S,i 


for  i  =  N1  +  2,  . ..,  1^  +  N  +  1  (8) 

n+l  n  ,  ^2^  /  n+l  n+l, 

ai  =  ai  +  ~E ^  (vi  ~  vi-i} 

for  i  =  N  +  2,  ...,  Nx  +  N2  +  2  (9) 

where  FQn.  indicates  that  the  displacement  un  is  used  in  the  force-displacement 

relation  Fc  =  -Ku  of  the  shear  support  unless  the  plastic  limit  F  controls, 
b  o 

From  the  boundary  conditions 

a?  -  -  P(t  )  (10) 

1  n 

where  P(t)  is  the  pressure  loading  at  x  =  0;  at  the  free  surface,  x  =  (H^  +  H^) 

CrN1+N2  +  2  =  °  (11) 

The  velocity  at  the  soil-structure  interface  is  handled  separately  by 

n+l  n  .  At  n  _n, .  2  „  n 

v.  =  v.  +  —  (a..,  -  o l„  . 

1  1  m  l+l  l  2mH  A  S,i 


for  i  =  +  1 


where  m  =  (p,Ax,  +  p„Ax„)/2. 


2.3  Parametric  Studies 


The  primary  set  of  parameters  for  this  model  are 
Geometry:  H^, 

Soil  properties:  p^,  M^, 

Structural  properties:  p M2 
Shear  support:  K,  Fq 


5 


Finite  difference:  =  40,  ^  = 

Artificial  viscosity:  a^/2  = 

ai+l/2  =  for  *  =  1>  2»  • • • » 

3nd  ai+1/2  =  0  for  i  ■  Hx+  1,  ....  Nx+N2+1 

The  remaining  parameters  are  varied  as  follows 

U  -  1,  2,  3,  4 
o 

k  =  0,  1,  10  (kips/in^) 

Elastic  Shear  Support  (F  =  °°) .  Stress  histories  at  x  -  20.3  feet  for 

cases  U  =  1,  2,  3,  4  with  k  =  0  (no  shear  support)  are  shown  in  Fig.  5, 
o 

and  the  corresponding  velocity  histories  at  the  interface,  x  =  20  feet, 
are  shown  in  Fig.  6.  Stress  histories  at  x  =  20.3  feet  for  U  =  1  with 
k  =  0,  1,  10  are  shown  in  Fig.  7,  and  the  corresponding  velocities  at 
x  =  20  feet  are  shown  in  Fig.  8.  Stress  histories  at  x  -  20.3  feet  for 
Uq  =  4  with  k  =  0,  1,  10  are  shown  in  Fig.  9,  and  the  corresponding  velocities 
at  the  interface  are  shown  in  Fig.  10. 

Plastic  Shear  Support.  Results  corresponding  to  plastic  behavior  of  the 
shear  support  are  presented  for  various  values  of  f  so  that  a  ductility  of  10 
is  achieved  (iteratively).  The  ductility  is  defined  as  the  ratio  of  maximum 
displacement  of  the  plastic  shear  support  divided  by  the  displacement  at  the 
elastic  limit  of  the  support  (fQ/k).  Stress  and  velocity  histories  are 
shown  in  Figs.  11  through  14  for  the  following  cases. 


Fig.  11: 

U  = 
o 

1, 

k  =  1,  f  =  0.246 

Fig.  12: 

U  = 
o 

1, 

k  =  10,  f  =  0.885 
0 

Fig.  13: 

U  = 
o 

4, 

k  =  1,  f  =  0.122 

0 

Fig.  14: 

U  = 
o 

4, 

k  =  10,  f  =  0.505 

10 
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.  UNCOUPLING  METHOD  -  ONE  DIMENSIONAL  STUDIES 


The  uncoupling  procedure  to  be  studied  here  is  a  special  case  of  a  more 
general  uncoupling  approach  described  in  Appendix  A.  In  the  present 
approximation,  the  calculation  is  separated  in  two  parts: 

(1)  Generation  of  free-field  (no  structure)  velocities  and  tractions 
on  the  surface  located  where  the  external  surface  of  the  structure 
should  be. 

(2)  Structural  analysis  using  the  free-field  data  from  (1)  in 
conjunction  with  the  plane  wave  approximation,  where  the  unload¬ 
reloading  impedance  (p1.C  )  is  used  for  the  soil.  Let  a  and 

XU  F 

Vp  denote  the  free-field  stresses  and  velociites,  respectively. 

The  response  of  the  structure  is  analyzed  with  the  boundary  condition 

°B  *  °F  +  PlCU(VB  ‘  V  <18) 

being  used  at  the  surface  of  the  structure,  where  a  and  v 

D  B 

denote  the  structural  boundary  stresses  and  velocities,  respectively. 
3.1  Free-Field  Data,  Op,  vf 

Free-field <po  structure)  data  were  calculated  by  the  finite  difference 

method  of  Eqs.  5,  6,  7  with  a  soil  depth  which  is  sufficiently  long 

(~40  feet)  that  no  reflections  from  below  could  affect  the  solution  at 

x  =  20  feet.  To  reduce  errors  associated  with  numerical  approximations,  the 

finite  difference  parameters  were  made  identical  to  those  used  in  the  soil- 

structure  interaction  calculations  of  the  preceding  section.  Hence, 

Axp  *  Hp/(Np  +  1/2)  with  H^  =  20  feet,  N^  =  40;  a  =  0.04  for  all  points. 

Stresses  and  velocities  at  x  =  20  feet  for  U  =  1  and  U  =  4  are  shown  in 

o  o 

Figs.  15  and  16.  Stresses  at  x  =  20  feet  were  obtained  by  averaging  the  two 
stresses  on  either  side  of  the  particle  velocity  at  this  depth. 


3.2  Structural  Analysis 

The  response  of  the  structure  was  analyzed  by  the  finite  difference 
method  of  Eqs.  7,  8,  9,  11,  12,  with  the  soil  mass  ignored  in  Eq.  12; 
the  structural  boundary  stress  in  Eq.  12  was  based  on  Eq.  18 

°i  =  °F  +  PlCU(vi  “  VF)  for  1  =  Ni  +  1  (19) 

The  accuracy  of  this  uncoupling  approximation  can  be  evaluated  by  studying 
the  calculated  velocity  time  history  at  the  boundary,  since  the  response  of 
the  structure  is  completely  determined  by  this  quantity. 

Figure  17  shows  the  velocity  time  history  for  the  case,  U  =1, 

o 

k  =  0,  (no  shear  support)  with  the  corresponding  soil-structure  interaction 

results  from  Fig.  8.  This  is  simply  the  plane  wave  approximation  for  the 

elastic  case  and,  as  expected,  the  comparison  is  excellent.  Similarly,  the 

comparisons  in  Fig.  18  for  U  ■  1,  k  ■  1,  f  =  °°,  and  in  Fig.  19  for 

U  =  1,  k  =  10,  f  =  °°,  are  both  excellent, 
o  o 

Figure  20  shows  the  comparison  for  Uq  =  4,  k  =  0,  fQ  =  00  with  the 

corresponding  soil-structure  interaction  results  from  Fig.  10.  Similarly 

the  comparison  for  U  •  4,  k  *  1,  f  *  00  is  shown  in  Fig.  21;  the  comparison 

for  U  =  4,  k  =  10,  f  =  00  is  shown  in  Fig.  22.  In  these  cases,  the  soil 
o  o 

is  hysteretic  but  the  structure  is  elastic. 

Figure  23  shows  the  comparison  for  U  “  1,  k  «  1,  f  =  0.246  with  the 

results  from  Fig.  12.  The  comparison  for  U  =  1,  k  =  10,  f  =  0.885  is  shown 

o  o 

in  Fig.  24.  In  these  cases  the  structure  responds  plastically  but  the  soil 
is  elastic. 

Figure  25  shows  the  comparison  for  Uq  =  4,  k  =  1,  f  =  0.122  with  the 
results  from  Fig.  14.  The  comparison  for  Uq  =  4,  k  =  10,  fQ  =  0.505  is  shown 
in  Fig.  26.  In  these  cases  the  soil  is  hysteretic  and  the  structure  behaves 
plastically. 
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3.3  Single-Degree-of-Freedom  Results 


The  finite  difference  calculations  in  Section  2  and  in  this  section 

so  far  were  done  with  five  stress  points  in  the  structure  (Nj  =  5).  However, 

it  is  interesting  to  note  that  waves  propagate  relatively  fast  in  the  structure 

so  that  the  velocity  time  histories  in  Figs.  17  through  26  can  be  reproduced 

quite  well  by  representing  the  structure  as  a  single-degree-of- freedom 

particle  with  velocity  v  and  mass  per  unit  area  m  =  p2^2  •  The  equation 

of  motion,  including  the  uncoupling  approximation,  then  becomes 

F 

mv  =  -[0F  +  -  v)  ]  +  —  ^ 20) 

Figures  27  through  30  show  selected  comparisons  of  velocity  time  histories 
obtained  using  Eq.  20,  with  corresponding  case  results  for  ^  =  5.  The 
differences  are  quite  small  even  in  the  nonlinear  problems. 


3.4  Variation  of  Soil  Impedance 


All  results  so  far  have  been  based  on  a  constant  value  of  for 

the  soil  impedance  in  the  uncoupling  approximation.  This  choice  can  be 
seen  to  be  a  good  one  for  the  present  studies  by  examining  Fig.  31  where 
the  results  of  using  the  loading  impedance  p  C  are  compared  with  the 

J.  Li 

corresponding  case  results  using  p^C^  an^  the  interaction  results  from 
Section  2.  The  very  early  time  results  using  p  C  compare  better  with 

1  L 

the  interaction  results  but  once  unloading  begins,  the  difference  becomes 
noticeable.  Still,  it  is  interesting  to  note  that  even  though  the  soil 
impedance  was  changed  by  a  factor  of  two,  the  later-time  error  is  only 
about  thirty  per  cent.  Thus,  the  choice  of  the  soil  impedance  seems  to  be 
less  critical  than  one  might  expect. 
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Other  ways  of  choosing  the  soil  impedance  in  the  uncoupling  approximation 

may  be  considered.  For  example,  it  might  seem  that  the  impedance  could  be 

made  to  depend  upon  the  state  of  t’  soil  by  first  choosing  p  C  when  a 

1  L  F 

initially  loads  the  structure,  then  switching  to  p  C  once  0  indicates 

X  U  F 

unload-reloading  behavior  is  taking  place  in  the  soil.  This  has  been  tried 
and  the  results  are  not  substantially  altered  by  this  procedure  except  for  the 
introduction  of  some  high-frequency  numerical  oscillations.  These  oscillations 
are  the  result  of  a  change  in  acceleration  whose  magnitude  depends  on 
PjCC  -  C^) ,  the  current  magnitude  of  v,  and  the  time  increment  At.  Such 
numerically  activated  accelerations  are  quite  undesirable  and  should  be 
avoided.  Smooth  variations  of  the  impedance  may  be  used  but  since  the  initial 
loading  phase  is  so  brief  for  explosive  loadings,  the  use  of  p^C^  at  all 
times  seems  to  be  a  reasonable  approximation,  at  least  until  further  study 
is  made.  If  the  soil's  unload-reloading  behavior  is  known,  then  a  smoothly 
varying  impedance  may  be  considered.  (However,  in  many  situations,  p^C^  is 
only  known  roughly  and  a  constant  value  may  be  acceptable.)  Additional  study 
may  be  needed  in  this  area,  see  Ref.  [5]. 
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.  MULTI-DIMENSIONAL  CONSIDERATIONS 


The  multi-dimensional  generalization  of  the  uncoupling  approximation, 

Eq.  18  (see  also  Appendix  A,  Eq,  A-6) ,  can  be  written  as 

!b  =  £f  +  P£  •  (~B  '  ~F}  (21) 

where  t„  and  v_  are  free-field  traction  and  velocity  vectors  at  the 

r  ~r 

surface  where  the  structure  should  be,  and  tg  and  Vg  are  traction  and  velocity 
vectors  at  the  boundary  of  the  structure.  From  wave  propagation  considerations, 
the  form  of  the  tensor  c  can  be  expressed  as 


c  =  + 


c  0 
s 


0  0 


where  c  is  given  in  a  local  coordinate  system  whose  first  coordinate 
corresponds  to  the  outward  normal  of  the  surface  of  the  structure,  and 
the  sign  is  chosen  to  produce  outward  wave  propagation  (radiation  damping); 
0^  and  cg  are  the  speeds  of  dilatational  and  shear  waves,  respectively, 
in  the  soil.  As  mentioned  in  the  previous  section,  a  smooth  variation  in 
the  unload- reloading  behavior  can  be  introduced  into  c,  if  necessary. 

The  geometry  of  a  particular  situation  can  affect  the  way  in  which 
the  uncoupling  approximation  is  used.  In  Fig.  1,  for  example,  waves  will 
be  reflected  back  and  forth  between  the  roof  and  the  ground  surface.  This 
effect  would  have  been  noticeable  in  the  results  of  the  previous  section 
if  the  calculations  had  been  carried  out  much  further  in  time.  This  effect, 
while  noticeable  in  the  one-dimensional  case,  is  not  likely  to  be  important 
in  two-  and  three-dimensional  problems  where  the  depth  of  the  roof  is  on 
the  order  of  one  or  more  times  the  horizontal  span  of  the  structure.  The 
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importance  of  such  geometrical  effects  must  be  considered  carefully  in  a 
particular  analysis  since  the  free-field  data  do  not  include  these  effects. 


A  refinement  of  the  uncoupling  approximation  can  be  used  to  include 
geometrical  effects,  if  necessary.  First,  instead  of  a  free-field  calculation, 
a  soil-structure  interaction  calculation  is  made  with  a  coarse  model  of 
the  structure,  and  tractions  and  velocities  are  recorded  at  points  on  the 
exterior  of  the  structure.  These  tractions  and  velocities  are  then  used 
in  place  of  the  free-field  data  in  Eq.  21  with  a  suitably  refined  model  of 
the  structure.  It  is  also  possible  to  record  tractions  and  velocities 
on  the  surface  of  an  imaginary  convex  surface  outside  the  structure  and 
then  use  a  refined  model  of  the  soil  and  structure  inside  this  surface. 

The  major  drawback  to  using  this  approach  is  that  the  first  stage  calculation 
will  generally  be  more  complicated  and  hence  more  costly  than  a  free-field 
calculation.  Detailed  studies  involving  multi-dimensional  situations  are 
in  progress  and  will  be  reported  separately. 
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5.  SUMMARY  AND  CONCLUSIONS 

An  uncoupling  approximation  for  analyzing  explosively  loaded  structures 
embedded  in  hysteretic  media  has  been  presented.  The  method,  based  on 
wave  propagation  considerations,  can  be  viewed  as  the  plane  wave  approximation 
extended  to  nonlinear  media.  Its  accuracy  has  been  shown  to  be  quite 
good  for  a  series  of  one-dimensional  soil-sturcture  interaction  problems, 
even  when  the  structure  and  surrounding  soil  are  both  nonlinear. 

It  should  be  noted  that  both  free-field  traction  and  velocity  vectors 
must  be  known  or  calculated  in  order  to  apply  the  method.  In  addition, 
the  soil's  unload- reloading  behavior  is  required. 

Multi- dimensional  problems,  and  techniques  for  improving  the  method, 
were  also  described  briefly.  Further  development  and  experience  are  needed, 
but  it  appears  that  reasonably  accurate  results  can  be  obtained  by  judicious 
application  of  these  techniques. 
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APPENDIX 


INTERACTIVE  SCHEMES  FOR  NONLINEAR  MEDIA 

Consider  a  structure  S  partially  or  totally  embedded  in  a  nonlinear 

medium  M.  Let  B  be  the  surface  common  to  S  and  M.  Assume  that  free  field 

input  is  available  and  could  be  used  in  a  simpli ried  structure-medium 

calculation  in  which  only  a  coarse  approximation  of  the  structure  is  used. 

T 

This  solution  gives  trial  boundary  tractions  t  . 

~B 

Let  the  correct  boundary  tractions  for  the  actual  structure  be  denoted 

T 

by  tg  .  Due  to  the  error  tg  -  tg  in  the  tractions  in  the  simplified  calculation, 

T 

the  velocities  along  B  are  in  error  by  an  amount  v  -  v  as  determined  by 

B  B 

Eb  -  £»  ‘  !{Eb  -  Eb}  <a-« 

where  F  is  a  functional  which  describes  the  response  of  the  surrounding 
material  to  an  excitation  on  boundary  B.  The  proposed  uncoupling  procedure 
consists  of  finding  an  appropriate  approximation  to  functional  F  . 

For  a  general  nonlinear  structure,  finite  element  equations  will  be 
of  the  form 

<Kx)  =  P  (A-2) 

in  which  x  and  P  are  the  displacement  and  interface  force  vectors,  respectively, 

and  <j>  is  an  operator  which  represents  the  structural  characteristics.  If 
T 

tg  and  tg  are  first  Piola-Kirchhof f  tractions,  P  depends  linearly  on  the 
interface  tractions  tfi,  i.e., 

P  -  P(tJ)  -  P(tg  -  tB)  (A- 3) 

Then  (A-2)  may  be  written  as 

4>(x)  +  P[F(Vg  -  yB)  ]  =  P(tg)  (A— 4) 
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V.V- 


The  left  hand  side  of  (A-4)  represents  a  modified  structural  operator  in 

which  the  boundaries  of  the  actual  structure  are  supported  by  nonlinear 

T 

hysteretic  connections  to  a  base  which  moves  with  velocity  v_  .  The  right 
hand  side  represents  the  interface  loading  from  the  trial  calculation. 


In  physical  terms  this  means  that  the  traction  t  can  be  used  to  load 

T 

the  actual  structure.  However,  to  account  for  the  fact  that  t  corresponds 

~D 

to  a  slightly  different  structure,  the  boundary  supports  of  the  actual 

structure  are  modified  by  the  use  of  non-linear  supports.  These  supports 

T 

are  attached  to  a  base  with  the  approximate  motion  v_  to  the  structure  at 
each  of  its  nodes  or  velocity  points.  The  support  behavior  F  is  chosen  so 
as  to  represent  the  appropriate  loading/unloading/reloading  behavior  of  the 
medium  surrounding  the  structure. 


As  a  first  choice,  take 

F  -  pc  (A-5) 

so  that  (A-l)  becomes 

*B  •  5b  =  P£  •  (?B  ■  V  (A'6) 

in  which  c  is  a  tensor  formea  from  the  wave  speeds  in  the  medium.  This  is 
equivalent  to  the  plane  wave  approximation  used  for  linear  media  and  involves 
only  linear  dashpot  supports  on  the  structure.  Similarly,  F  can  be  chosen 
to  give  an  interactive  scheme  analogous  to  the  doubly  asymptotic  approximation. 
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AIR  OVERPRESSURE 


FIG  I.  A  SOIL- STRUCTURE  INTERACTION  PROBLEM 


FIG  2.  ONE  DIMENSIONAL  SOIL-STRUCTURE  INTERACTION  MODEL 
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FIG.  5  STRESS  HISTORIES  AT  20.3  FEET  FOR  SOIL 

UNLOADING  CASES  U.=  l,2,3,4  WITH  k  =  0,  f„=  00 
(NO  SHEAR  SUPPORT) 
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FIG.  6  VELOCITY  HISTORIES  AT  20  FEET  FOR  SOIL 

UNLOADING  CASES  U0  ■  1,2, 3,4  WITH  k«0,f.»00 
(NO  SHEAR  SUPPORT)  23 
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FIG.  9  STRESS  HISTORIES  AT  20.3  FEET  FOR  HYSTERETIC 
SOIL  IU«=4)  AND  ELASTIC  SHEAR  SUPPORTS  (k=0,l,IO) 
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FIG.  10  VELOCITY  HISTORIES  AT  20  FEET  FOR  HYSTERETIC 
SOIL  (U.*4)  AND  ELASTIC  SHEAR  SUPPORTS 
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FIG.  II  STRESS  HISTORIES  AT  20.3  FEET  FOR  ELASTIC 
SOIL  (Uo=l)  AND  PLASTIC  SHEAR  SUPPORTS. 
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FIG.  12  STRESS  HISTORIES  AT  20  FEET  FOR  ELASTIC 
SOIL  (U.*l)  AND  PLASTIC  SHEAR  SUPPORTS 


STRESS 


k  «  l 
f*  =0.122 


/  k.io 

/  fo«  0.506 


t,=  10 


*  -8.60 


TIME  MSEC 

FIG.  13  STRESS  HISTORIES  AT  20.3  FEET  FOR  HYSTERETIC 
SOIL  (U.«4)  ANO  PLASTIC  SHEAR  SUPPORTS 
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FIG.  14  VELOCITY  HISTORIES  AT  20  FEET  FOR  HYSTERETIC 
SOIL  (U«*4)  ANO  PLASTIC  SHEAR  SUPPORTS 
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FIG.  16  FREE  FIELD  VELOCITIES  AT  20  FEET  FOR  ELASTIC 
(U.*l)  AND  HYSTERETIC  ( U©  s  4)  SOILS. 
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COMPARISON  OF  RESULTS  FROM  UNCOUPLING 
APPROXIMATION  WITH  SOIL- STRUCTURE  INTERACTION 
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FIG.  17 

ELASTIC  (U.  =  I)  SOIL, 

NO  SHEAR  SUPPORT  (  k  =0) 


IS. 00  20.00 

TIME  MSEC 


2S.00 


FIG.  18 

ELASTIC  (  Uo  s  I)  SOIL, 
ELASTIC  tk  =  l)  SHEAR  SUPPORT 
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FIG.  19 

ELASTIC  (U*  =  I)  SOIL, 
ELASTIC  (k  =  IO)  SHEAR  SUPPORT 
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FIG.  20 
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FIG.  21 

HYSTERETIC  (U**4)  SOIL, 
ELASTIC  (k=l)  SHEAR  SUPPORT 
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FIG.  22 

HYSTERETIC  (U»=4)  SOIL, 
ELASTIC  (k *  10)  SHEAR  SUPPORT 
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FIG.  23 

ELASTIC  (Uo  =  l)  SOIL, 
PLASTIC  IK* I ,  to* 0.246) 
SHEAR  SUPPORT 
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FIG.  24 

ELASTIC  (U**l)  SOIL, 
PLASTIC  (k  =  10,  f0*  0.885) 
SHEAR  SUPPORT 
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FIG.  25 
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FIG.  26 
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FIG.  27 
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FIG.  28 
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FIG.  30 

HYSTERETIC  (U«  =  4)S0IL, 
PLASTIC  (k«*l,  foB  0.122) 
SHEAR  SUPPORT 
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COMPARISON  OF  UNCOUPLING  APPROXIMATION 
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FIG.  31  HYSTERETIC  (Uo  =  4)  SOIL,  NO  (k=0) 
SHEAR  SUPPORT 
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